Exact Non-Stationary Probabilities in the Asymmetric Exclusion Process on a Ring 
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By a geometrical treatment of the Bethe ansatz, we obtain an exact solution for the totally 
asymmetric exclusion process on a ring. We derive an explicit determinant expression for the 
non-stationary conditional probability Prob(x\, ...,xp\t\x1, ...,x%;0) of finding P particles on sites 



xi, ip at time t provided they are on sites x\, x° P at time t — 0. 



PACS numbers: PACS numbers: 05.40.+J, 02.50.-r, 82.20.- 



The one-dimensional asymmetric exclusion process 
(ASEP) has been intensively studied as one of the sim- 
plest examples of a system with stochastic dynamics 
and exclusion interaction 0],[|,g. In the steady state, 
many properties of the ASEP have been calculated ex- 
actly (see, e.g., [[|] • H , § . @] an d references therein). The 
structure of states out of stationarity is more compli- 
cated and description of dynamics is harder to obtain 
H , H . In general, evaluation of the conditional probabil- 
ity of finding P particles on lattice sites x\, xp at time 
t with initial occupation x\, ...,x P at time t = 0, needs 
solving a non-trivial master equation. The solution of 
the master equation for the ASEP with a finite num- 
ber of particles on an infinite lattice was obtained in |l(| 
where an explicit determinant expression for the proba- 
bility Prob{x\, xp; t\x\, x P ; 0) was derived. At the 
same time, a non-stationary solution for the finite density 
of particles is desirable for considering the evolution of 
a system from arbitrary initial conditions to the steady 
state. 

It is the aim of this Letter to present an exact expres- 
sion for the conditional probability of finding P particles 
on an finite ring of L sites after the evolution with the 
ASEP dynamics during finite time t. In the theory of 
exactly soluble models, the determinant solutions on fi- 
nite two-dimensional domains permitting a dynamical in- 
terpretation in space-time coordinates are known either 
for so-called free-fermion models pd[| , or for particular 
domain- wall boundary conditions ]12[,p^]. The ASEP 
does not belong to the free-fermion class and cannot be 
restricted to a domain wall geometry. However, we will 
see that a specific property of the pair interaction in the 
ASEP allows one to consider the model on the finite ring 
and express the solution in a form of determinant of a 
P x P matrix. 

Instead of explicitly solving the master equation on the 
ring, we use here a geometrical treatment of the Bethe 
ansatz which allows us to analyse a tangled system of 
allowed and forbidden trajectories of particles. To make 
the presentation more transparent, we consider a discrete 
space-time version of the totally ASEP. The continuous 



time limit then follows from the final expressions by the 
straightforward substitution of the Bernoulli distribution 
by its Poisson analogue. 

Consider the triangle lattice A obtained from the 
square lattice by adding a diagonal between the upper 
left corner and the lower right corner of each elementary 
square. The lattice is periodic in the horizontal direction 
with the period L. Let (x, t) be integer coordinates of 
a particle on A where the vertical time axis is directed 
down and the horizontal space axis is directed right. A 
trajectory of the particle is a sequence of connected ver- 
tical and diagonal bonds of A. Each diagonal bond cor- 
responds to one jump of the particle to its right for a 
unit of time and has a statistical weight z. The vertical 
bond corresponds to a stay at given site for the unit time 
interval and has a statistical weight y. The generating 
function of all possible free trajectories of one particle for 
time t is (z + y) . For trajectories starting at an initial 
position (x°,0) and ending at (x,i), we define 



B(N,t) 



JV 



(1) 



where N 



x° is the travelled distance. To provide 



B(N, t) with probabilistic meaning, we put z < 1, z + y = 
1. Then, B(N,t) is the probability to reach (x,t) from 
(x°,0) for t time steps provided z is the probability of 
one step in the right direction. In the continuum limit, 
z — > 0, t — ► oo, we have 



B(N,t) 
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(2) 



where t is rescaled continuous time tz — > t. For N < 0, 
we put B(N,t) = 0. 

The ASEP in discrete space-time can be defined as 
follows. Consider trajectories of P particles on the lattice 
A starting at points (x°, 0), (x° p , 0), < x? < x° < 
... < x P < L and ending at points (xi,t), (xp,t), < 
xi < X2 < ... < xp < L after an arbitrary number 
of rotations around the ring. The exclusion rules read: 
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(a) trajectories of particles do not intersect; (b) for an 
arbitrary elementary square of A, if two vertical bonds 
are occupied by adjacent trajectories, the weight of the 
left bond is changed from y to 1. The rule (a) is the 
usual condition of occupation of every site by at most one 
particle. The rule (b) implies that the moving particle 
stays at a given site with probability 1 if the target site 
is occupied by a standing particle. 

The generating function Gp(x\, xp; t\x{, x P ; 0) 
of the discrete ASEP is the sum over all trajectories of P 
particles allowed by the exclusion rule (a) and weighted 
according to the rule (b). Due to the exclusion rules, Gp 
is Prob(x\, xp; t\x\, x p ; 0). 

To formulate the main result of this Letter, consider 
the function JhJ 

F m (xl Xj \t) =£( fc t-7 1 ) k > x ^ (3) 

fc=0 ^ ' 

if integer m > 0, and 

F m {xlx s \t) = ( ~ m ) F (s? - k, Xj \t) (4) 

if integer m < 0. For m — 0, 

F (xlx j \t) = B(x j -x i ,t) (5) 

where B(N,t) is given by Eq. (Q) . Then, for integer 
L > 1 and < P < L the generating function 
G P (xi, ...,xp;t\xl, ...,x° P ;0) is 

Gp= J2 - J2 (-l)^'<^"^" 3l detM (6) 

ni— — oo np — — oo 

Elements of the P x P matrix M are 

M ij = F aij (ril,x j + n j L\t) (7) 

where 

Sij = (P - l)nj -y^k+j-i (8) 

In the continuous limit, B(N,t) is given by 
Eq.(||) and Gp coincides with the probability 
Prob(x 1 ,...,x P ;t\x^,...,x P ;0) of the totally ASEP 
in its standard formulation. 

The derivation of this result is based on a com- 
mon property of integrable models admitting of a two- 
dimensional graphic representation: interchanging of end 
points of two trajectories leads to their crossing. The 
idea of the Bethe-ansatz is to represent trajectories of 
interacting particles by a set of free trajectories given 
by Eq.(g) or Eq.(g). Then, using the one-to one corre- 
spondence between intersections and permutations, one 



can reduce enumerating all interacting trajectories to a 
proper choice of signs of permutations. 

We start with the case of two particles P = 2. Ac- 
cording to the Bethe ansatz, we try to represent the 
motion of particles by free trajectories from (x°,0) to 
(xi,t),i — 1,2. Consider an elementary square of A with 
space coordinates x of the left side and x + 1 of the right 
side. Assume that particles come for the first time to 
neighboring sites at a moment t when one trajectory 
reaches the site (x,t ) from (x?,0) and another reaches 
the site (x+l,t ) from (x°, 0). To ensure correct weights 
of the next steps of interacting particles after moment t , 
we have to exclude two possibilities from all continuations 
of trajectories (Fig. la): 

(i) for the first particle,the step from (x, t ) to (x+l, t + 
1) with weight z and then from (x + l, t +1) to (xi, t); for 
the second particle, the step from (x+l, t ) to (x+l, t +1) 
with weight y and then from (x + l,t + 1) to (x%, t). 

(ii) for the first particle,the step from (x, t ) to (x, t +1) 
with weight y — 1 = — z and then from (x, t + 1) to 
(x\,t); for the second particle, the step from (x+l,t ) to 
(x + l,t +1) with weight y and then from (x + 1, t +1) 
to (x%,t). 




FIG. 1: The interaction between two trajectories (see text). 

Case (i) is the forbidden step of the first particle toward 
the standing second particle. Case (ii) is a correction of 
the weight of the vertical step of the first particle which 
must be 1 instead of y according to the ASEP rule (b). 
The generating function of paths of the first particle in 
case (i) is a product of three factors B(x — x\,t )zB(x\ — 
x — 1, t — t — 1). Using symbolic notations W(a, x\z\x + 
1,6) for the generating function of trajectories passing 
points a, x, x + 1, b at moments 0, t , t + l,t and making 
a step with weight z between t and t +1, we can write 
the contribution from diagram (i) in the form 

Wi = W(x\,x\z\x + l,xi)W(x° 2 ,x + l\y\x + l,x 2 ) (9) 

The contribution from the diagram (ii) is 

W 2 = -W(x° 1 ,x\z\x,x 1 )W(x° 2 ,x + l\y\x + l,x 2 ) (10) 

Consider now the trajectories where the end points are 
interchanged (Fig lb). The contribution from these di- 
agrams is W(x\, x\z\x + 1, x 2 )W(x 2 , x + l\y\x + 1, xi) — 
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W(xi,x\z\x,x 2 )W(x 2l x + l\y\x + l,xi). We are going to 
take the diagrams in Fig. lb with opposite signs to can- 
cel W\ + W 2 . The left diagrams in Fig. la and Fig. lb 
are equivalent, however the right ones are different. To 
cancel all unwanted diagrams, we add to the diagrams 
in Fig. lb a set of auxiliary trajectories. Namely, add to 
trajectories of the second particle those starting in point 
x® — 1 and taken with minus sign. Also, we add to trajec- 
tories of the first particle a set of trajectories starting at 
the points shifted along the ring in negative directions: 
x\ — l,x\ — 2, x\ — 3,... . If a shift exceeds kL, where 
k > is integer, the trajectory wraps the cylinder A k 
times. Then, the contribution from diagram (i) in Fig. lb 
will be 



(11) 



where 



W£ = ^W{x\ -k,x-k\z\x + l -k,x 2 ) (12) 



fc=0 



and 



Wf = W(x%,x+I\y\x+I,xi)-W(x%-I,x\y\x,xi) (13) 
Correspondingly, for the diagram (ii) in Fig. lb, we have 
W 2 = W^W 2 (14) 

where 



at least once the point (a; — k, t ) or (x + 1 — k, t ) where 
it participates in the cancellation procedure. 

Each free trajectory from a to b making the vertical 
step at the collision site x can be decomposed into two 
parts W(a, x\l\x, b) + W(a, x\ — z\x, b). The second part 
is unwanted and is cancelled, the first one corresponds to 
trajectories which continue with true weights up to the 
next collision. As the second term in Eq.(]l7|) contains 
intersecting trajectories only, all of them will be cancelled 
eventually and only true allowed trajectories from the 
first term survive. 

The ASEP on the ring has several peculiarities. To fix 
them, let us note that two intersecting trajectories are 
non-equivalent: one of them belongs to the overtaking 
particle and we may call it " active" . On the contrary, 
the second particle can be called "passive". In the case 
of infinite lattice, the active and passive trajectories are 
ordered: for each pair i, i + 1, the trajectory of i-th par- 
ticle with respect to i + 1 particle is always active. On 
the ring, each of two trajectories can be active or passive 
independently on initial conditions. Moreover, one tra- 
jectory can intersect another m times if the numbers of 
rotations differ by m for two particles. 

Assume, that the trajectory of given particle has m ac- 
tive intersections. It means that it participates m times 
in the cancellation procedure and its starting point is 
shifted m times to arbitrary distances in the negative 
direction of the ring. As a result, the auxiliary set asso- 
ciated with the free trajectory between points x® and Xj 
becomes 



W. 2 + = -^2w(x° 1 -k,x-k\z\x-k,x 2 ) (15) 

and W 2 — Wi . Taking into account that generating 
functions of trajectories from — k,0), i = 1,2, to 
(x — k.t ) are equal for all k due to translation invariance, 
one can check the identity 



Wi + W 2 - Wi - W 2 = 



(16) 



comparing all positive and negative terms. 

Consider first the evaluation of generating function 
G 2 (xi,x 2 ]t\x\,x 2 ]Q) in the case L ^> t and L ^> x 2 > 



3> which is eq 



uivalent to the ASEP on an infinite 



lattice solved in |fuj |. In this case, 

G 2 = B{xi - x\, t)B(x 2 -x%,t) -{B{ Xl -x° 2 ,t) 
-B{x x -x° 2 + l, *)) J2Zo B ^ -4 + k, t) 

(17) 

Indeed, the first term in Eq.(|l~7|) generates all possible 
free trajectories from initial to end points. When one par- 
ticle approaches another, the second term produces tra- 
jectories cancelling unwanted terms. On the other hand, 
the order of starting and ending points in the second 
term is interchanged. Therefore, each trajectory from 
the second term starting at x\ — k or x\ — k approaches 



B{x 3 ~xlt)^Y. 



k=0 



k + m — 1 

771—1 



B(xj-x^+k,t) (18) 



because the shift by k positions for m attempts can be 
done in (k + m — l)!/(m — l)\k\ ways. The above can be 
expressed in an operator form 



1 



(i - a*)' 



-B(x 3 -xlt) (19) 



where the operator a, shifts x® by one step in the nega- 
tive direction. Similarly,for trajectories having m passive 
intersections we get 



k=0 



Bfa -x° t +k, t) (20) 



because the right-hand side is the result of action of the 
operator (1 — ai) m . Note that Eq.(p^[) coincides with 
Eq.(||) and Eq.(|20|) with Eq.(Q) where the index m can 
be called " activity" . 

To find the generating function of two particles on the 
ring, we map the trajectories wrapping the cylinder A on 
an infinite plane introducing coordinates x + nL, n inte- 
ger, for equivalent points. We call trajectories of two par- 
ticles compatible if there is at least one possibility to draw 
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them without intersections. Given starting points and 
x°, at t — 0, the pairs of compatible trajectories corre- 
spond to end points x\, Xi\ X2, x\ +L;...;xi + nL, X2 + nL; 
X2 + nL, xi + (n + 1)L;... at time t. The index of activ- 
ity of compatible trajectories is 0. A trajectory ending 
at x\ + nL may interact with the trajectories ending at 
X2 + {n — l)L or X2 + nL. Following the Bethe-ansatz 
prescription, we should add to the set of compatible tra- 
jectories two sets of intersecting trajectories taken with 
minus sign: the set which is obtained by permuting end 
points X2 + (n — 1)L and x\ + nL and the second one 
obtained by permuting x\ + nL and X2 + nL. Each new 
trajectory, in its turn, interacts with neighbouring tra- 
jectories and we should permute their end points again. 
Continuing this procedure, we obtain all possible pairs of 
trajectories of the first particle with end point x\ + n±L, 
the second one with end point X2 + n2L, or vice versa, 
for arbitrary integer n\ and «2 • The permutations of end 
points can be expressed by the determinant as can be 
seen m Eq.(g). The sign of a pair is defined by the num- 
ber of permutations needed to obtain given pair from an 
compatible one. The index of activity of each trajectory 
is defined by the number of overtakes to the moment t. 
Evaluation of the number of permutations gives the pre- 
factor in Eq.(^). The number of overtakes is given by 
Eq.®. 

If the number of particles P > 3, the elementary 
squares shown in Fig.l may occur several times in one 
horizontal strip of A. If squares filled by interacting 
trajectories are separated one from another by a gap 
of empty sites, the above arguments can be applied to 
each pair of interacting trajectories separately. The cru- 
cial case for the Bethe ansatz is a situation when the 
elementary squares are nearest neighbors. The specific 
property of the totally ASEP is that, in each pair of in- 
teracting trajectories, the right trajectory remains free 
and interacts with the next trajectory independently on 
its left neighbors. Therefore, we can analyse the inter- 
action between particles considering successively elemen- 
tary squares in each row from left to right starting from 
an arbitrary empty square and then from the top to bot- 
tom of the lattice until all unwanted trajectories on A will 
be removed. As above, all trajectories ending at points 
Xi +mL,i = 1, P must be involved in consideration. 



Then, Eq.(^) is a straightforward generalisation of the 
case P = 2, and Eq.(^)gives the index of activity for an 
arbitrary number of intersecting trajectories. A new el- 
ement in the many-particle case is that each trajectory 
may get m active intersections and n passive ones. The 
operator form of Eq. ( [l9| ) and Eq. ( p0| ) shows that the re- 
sulting activity in this case is m — n. 

Because of conditions B(N,t) = for N < 0, the in- 
finite summations in Eq.(j^) and Eq.(|6[) are actually fi- 
nite for finite t. An advantage of the discrete formula- 
tion of the ASEP is a possibility to illustrate each step 
of derivation by simple examples. For instance, in the 
case P = 2, L = 3, t = 6, xi = xi = l,^ — x i = 2, 
the generating function is G2 — z 12 + 20z 9 y 2 + 30z 9 + 
z 6 y 6 + 30zV + 90zV + 20zV + 20z 3 y 6 + 30z 3 y 5 + y 6 . 
All 243 allowed configurations enter this expression with 
proper weights. In the general case, the obtained ex- 
pressions open a prospect for detailed investigations of 
non-stationary solutions for various initial conditions. 
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